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Abstract 

We consider the physical combinatorics of critical lattice models and their associated conformal 
field theories arising in the continuum scaling limit. As examples, we consider A-type unitary 
minimal models and the level- 1 s£(2) Wess-Zumino-Witten (WZW) model. The Hamiltonian of 
the WZW model is the U q (s£(2)) invariant XXX quantum spin chain. For simplicity, we consider 
these theories only in their vacuum sectors on the strip. Combmatorially, fermionic particles 
are introduced as certain features of RSOS paths. They are composites of dual-particles and 
exhibit the properties of quasiparticles. The particles and dual-particles are identified, through 
a conjectured energy preserving bijection, with patterns of zeros of the eigenvalues of the fused 
transfer matrices in their analyticity strips. The associated (m, n) systems arise as geometric 
packing constraints on the particles. The analyticity encoded in the patterns of zeros is the key 
to the analytic calculation of the excitation energies through the Thermodynamic Bethe Ansatz 
(TBA). As a by-product of our study, in the case of the WZW or XXX model, we find a relation 
between the location of the Bethe root strings and the location of the transfer matrix 2-strings. 
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1 Introduction 



The phrase physical combinatorics [H [2] was termed in 2000 to refer to a new body of results ap- 
pearing at the intersection of solvable lattice models [3], representation theory of Conformal Field 
Theory (CFT) and combinatorics. The term quasiparticle takes somewhat different meanings ac- 
cording to the field, but we use (fermionic) quasiparticles in the sense described by McCoy and 
his collaborators [1] in the early nineties. Quasiparticles have properties usually associated with 
particles but may not correspond to real physical particles. More specifically, physical combina- 
torics is the combinatorial realization of finitized conformal characters El (or more generally 
conformal partition functions built from such characters) as generating functions of diagrammatic 
objects weighted by an energy statistic E€N for system size N 

x2°(ff) = T c/24+A £ Q E ^ ect \ XA(q) = Jim X f\ q ) (1.1) 

' — ' N— >oo 

object 

where c is the central charge. In the thermodynamic limit N — > oo, this leads to fermionic 
expressions for the conformal characters xa(q) as positive term series in the modular nome q = e nlT 
where r is the modular parameter. In the limit q — > 1, the finitized character Xa (?) counts the 
number of states for system of size N. The diagrammatic objects can be strings, lattice paths, 
rigged configurations, Young tableaux or other suitable combinatorial construct. Ultimately, these 
combinatorial objects should manifest the nature of the constituent fermionic quasiparticles. 

Historically, physical combinatorics first appeared in the context of the string hypothesis [8\ of 
Bethe's solution [9] of the Heisenberg (spin-^ XXX) quantum spin chain. For each allowed eigen- 
value of the commuting family T(u) of finite-size transfer matrices, Takahashi's string hypothesis 
posits the allowed patterns of zeros in the complex plane of the spectral parameter u. More cor- 
rectly, the hypothesis specifies the zero patterns of the eigenvalues of Baxter's auxiliary matrices 
Q(u) but the eigenvalues T(u) are easily reconstructed from the eigenvalues Q{u) through the 
Bethe ansatz [9], or equivalently, Baxter's T-Q relation |10j . In the intervening years, there have 
many works [HI [121 d3l U31 USl UBl UZl USl USl EQl EH E21 E31 developing the enumeration and study 
of Bethe states in terms of strings, rigged configurations and Young tableaux. Most notably, in 
the case of the XXZ spin chain with periodic boundary conditions, there is still some debate about 
whether the set of Bethe states given by the string hypothesis is complete \24\ [25| 126] . 

A major advance in the methods of physical combinatorics came with the advent of Bax- 
ter's one-dimensional configurational sums X^ N \q) which derive from off-critical Corner Transfer 
Matrices (CTMs) [3J. The one-dimensional configurational sums are generating functions for one- 
dimensional lattice paths weighted by an energy statistic. In the prototypical example of the 
minimal models [27J with central charges 

C = 1 "I(TTT) < L2 > 

the conformal characters exactly coincide [5] with the N — > oo limit of the one-dimensional con- 
figurational sums of the Andrews-Baxter- Forrester (ABF) Restricted-Solid-On-Solid (RSOS) mod- 
els [2H] in the off-critical regime III. This correspondence principle further extends [29, El E] to 
finitized conformal characters 

X {N \<1) = q- c/2A X^(q) . (1.3) 

q=q 
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In applying physical combinatorics, some form of this correspondence principle is invariably invoked 
to introduce RSOS lattice paths through the one-dimensional configurational sums. However, 
substituting q with q is surprising and even logically inconsistent. Specifically, the modular nome 
q in the unitized characters is related to the aspect ratio and conformal geometry, and only makes 
sense at criticality, whereas the elliptic nome q in the one-dimensional configurational sums measures 
the departure from criticality in off-critical models. Reinterpreting q as q does not substitute 
for understanding how the structure of RSOS paths actually arises in integrable critical lattice 
models. While the approach based on CTMs and the correspondence principle has proved very 
successful [HE], it is natural to ask whether this approach to physical combinatorics really is 
physical. This approach also has the inherent limitation that it can only be applied with fixed 
and not periodic boundary conditions. There have been many subsequent works studying the 
eigenvalues and physical combinatorics of the RSOS models both with periodic [301 ED E2] and 
fixed [331 EH E3 EH EZ] boundary conditions but perhaps the RSOS lattice path approach to the 
physical combinatorics of the RSOS is best typified by the work of Warnaar [38} 139] . For some 
recent reviews see |40| . 

The aim of this paper is twofold — first, to present a conjectured energy-preserving bijection 
between RSOS paths and string patterns of the transfer matrices at the critical point and, second, 
to give complete derivations of the conformal energies in the vacuum sectors of the RSOS minimal 
and XXX models by solving the Thermodynamic Bethe Ansatz (TBA) [41, 32J. Since the solution of 
the TBA equations relies only on the string patterns, it is independent of the conjectured bijection 
with RSOS paths. Following initial work in [42] . the approach we adopt in this paper removes the 
logical inconsistency of invoking properties of off-critical models to study the physical combinatorics 
of critical RSOS models. It points the way to a stronger unification of the approaches based on 
string patterns and RSOS lattice paths. Moreover, our methods can be applied to larger classes 
of boundary conditions, including models with periodic boundaries which are not tractable by 
previous RSOS lattice path methods. 

In this paper we consider the ABF RSOS Al models and the s£(2)i WZW lattice model in 
their vacuum sectors. The WZW model is a six- vertex model with rational Boltzmann weights 
which reduces, in the Hamiltonian limit, to the C7 g (s£(2))-invariant XXX quantum spin chain. In 
the sequel, we will always refer to this case as the XXX model. The layout of the paper is as 
follows. Quasiparticles are introduced combinatorially by decomposing one-dimensional RSOS lat- 
tice paths into particles decorated by dual-particles in Section 2. It is shown that the geometric 
packing constraints of these decorated particles encode the (m, n) systems. Moreover, it is argued 
that the RSOS lattice paths actually encode string patterns through the locations of the quasipar- 
ticles (particles and dual-particles). This is argued on the basis of a conjectured energy-preserving 
bijection between the RSOS paths and the string patterns. The energy statistics are associated 
with the string patterns through unitized characters (involving integer string quantum numbers 
and the Cartan matrix) and with RSOS paths through one-dimensional configurational sums (hav- 
ing their origins in off-critical Corner Transfer Matrix calculations). In Section 3, we consider the 
example of the vacuum sector of the minimal RSOS lattice models. Using the commuting family 
of fused transfer matrices, we show how the eigenvalues are completely classified by string patterns 
consisting of 1- and 2-strings in L — 2 analyticity strips. The allowed numbers of 1- and 2-strings in 
these strips (particle content) is subject to the (m, n) system. The relative positions of the 1- and 
2- strings in each strip (integer string quantum numbers) uniquely specifies the analyticity input 
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required to analytically solve the TBA equations. The detailed solution of the TBA presented in 
this section finally leads to the general conformal energy expression (j3.136H including the central 
charge term. This is the first complete derivation of this general result from the TBA. In Section 4, 
we consider the vacuum sector of the XXX model. Although this model has not previously been 
analysed using RSOS paths, we show that this model is amenable to our new approach. Considering 
the commuting family of fused transfer matrices again shows that the eigenvalues are completely 
classified by string patterns consisting of 1- and 2-strings in L = \_N/2\ + 1 analyticity strips. In 
this case the fusion hierarchy is infinite and does not truncate. This is the reason that L grows 
with the system size N. The allowed numbers of 1- and 2-strings in these strips (particle content) 
is again subject to the (m, n) system given in Section 2. The relative positions of the 1- and 2- 
strings in each strip (integer string quantum numbers) uniquely specifies the analyticity input re- 
quired to analytically solve the TBA equations. The detailed solution of the TBA finally leads to 
the general conformal energy expression (|4.7ip including the central charge term. Again this is the 
first complete derivation of this general result from the TBA. We conclude with a discussion in 
Section 5. Details on dilogarithm identities [32} |4"31 EH] are relegated to an Appendix. 

2 Combinatorial Quasiparticles 

In this section we introduce quasiparticles [I] of two kinds in a purely combinatorial manner and 
refer to them as particles and dual-particles even though they are all in fact quasiparticles. The 
connection with various lattice models is established in the next section. 

Combinatorially, quasiparticles are introduced via features of Restricted Solid-On-Solid (RSOS) 
paths. In principle, the paths can be arbitrary paths on Dynkin diagrams of A-D-E-T or more gen- 
eral type, but here we only consider Al (linear) and Tl (tadpole) diagrams with L nodes (Figure [T]). 
We denote by Tl the tadpole diagram with a single loop at node j = 1 and by T' L the tadpole 
diagram with a single loop at node j = L. On the Dynkin diagram Al with L fixed, an TV-step path 
a = {ao, ai, a2, • • • , a^} satisfies the RSOS constraint aj_i — aj = ±1 with aj-i G {1, 2, . . . , L} for 
j = 1, 2, . . . , N. On the tadpole diagram Tl paths with steps along the baseline with Oj-i = a* = 1 
are also allowed. In this case, we take L = \_N/2\ + 1 corresponding to the maximum height oc- 
curring in an A-step path. We stress that this means, in the XXX case, L grows with the system 
size N. We will only consider the vacuum sector with boundary conditions given by ao = ajv = 1) 
dN+l = 2. 



A L minimal: G = A L I 1 1 1 1 1 1 G* = A L -2 I 1 1 1 1 1 1 

1 2 3 ••• L 1 2 3 L-2 

s£(2) 1 XXX: G = T L Q— I 1 1 1 1 1 G* = T' L _ X \ 1 1 1 1 1— Q 

12 3 ••• L 12 3 ••• L-1 

Figure 1: Dynkin diagrams G = Al,Tl and their respective duals G* = Al-2^T' l _ x . The particles 
live on G and the dual-particles live on G*. In the Al minimal case L is fixed. In the XXX case 
L = LA/2J + 1 grows with the system size N. 
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2.1 Pure particles 



On Al, pure particles |38|, [39] of type a correspond to upright two-dimensional pyramids (trian- 
gles) of height a and baseline width 2a. We use the name pyramid to emphasize that there is a 
distinguished base (baseline). There are thus L — 1 types of pyramid particles. In addition, on 
Tl, we allow tower particles consisting of a single step along the baseline with aj-\ = aj = 1. A 
typical path can contain many particles. Let n a with a = 1, 2, . . . , L — 1 be the number of pyramid 
particles of type a and let no denote the number of tower particles. Since there can only be or 1 
particle of a given type at a given j, there is an exclusion principle and the particles are fermionic 
in nature. 

n 4 re 3 n 2 n\ \n \ 
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1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 

Figure 2: Illustration of the different types of pure particles that can occur in an iV-step path. There 
are L — 1 types of pyramid particles of different heights. Their positions (peaks) are indicated by 
a solid coloured dot according to their type: • • • ■ • • A pyramid particle of type a is 

decorated by a(a — 1) dual-particles as shown. On Al, there are L — 2 types of dual-particles with 
different relative heights above the particle baseline. On Tl, there is an additional dual-particle 
giving a total of L — 1 different types of dual-particles. The dual-particles are also indicated by solid 
coloured dots according to their type: • • • • • • • In addition, on Tl, there are tower particles 

corresponding to single steps along the baseline with aj-i = aj = 1. These particles which are 
designated to be of type a = are decorated by a complete tower of L — 1 = [N/2\ dual-particles. 



2.2 Dual particles 

On Al, a particle of type a is decorated by a(a — 1) dual-particles as indicated in Figure 2. In 
this sense particles are composites of dual-particles. On Al and Tl, the type of the dual-particle 
is fixed by the relative height above the particle baseline. Specifically, a dual-particle of type a 
appears at a height a above the baseline of the particle which contains it. On Al, there are L — 2 
types of dual-particles whereas, on Tl, there are L — 1 different types. On Tl, the tower particles 
are decorated by a complete tower of L — 1 = L-^/^J dual-particles. A typical path can contain 
many dual-particles. Let m a with a = 1, 2, . . . , L — 1 be the number of dual-particles of type a. 
Since there can only be or 1 dual-particle of a given type at a given lattice point j, there is an 
exclusion principle and the dual-particles are also fermionic in nature. 

2.3 Geometric packing and (m, n) systems 

The (m, n) systems [61 [7] describe the possible particle contents. In our combinatorial framework, 
they arise from the geometric packing constraints of particles and dual-particles along an A^-step 



6 



(2.2) 



path. The adjacency matrices that appear in the (m, n) systems differ from the adjacency matrices 
of G = Al , Tl for the original paths and are in fact those of the dual graphs G* . Let A denote the 
adjacency matrix of G* = Al-2 or T' L _ X so that the corresponding Cartan matrix is 

C = 2I-A (2.1) 

Al: The geometric packing constraints to accommodate n a particles of type a and m a dual- 
particles of type a along a path of N steps on G = Al are 

N = 2ni + 4n 2 + 6n 3 + 8n 4 H h 2(L - l)n L -i 

mi = 2n 2 + 4n 3 + 6n 4 + • • • + 2(L — 2)ul-\ 

m 2 = 2n 3 + 4n 4 H h 2(L - 3)n L -i 

m 3 = 2n 4 H \-2(L- 4)n L _i 

rni-2 = 2n L -i 
Rearranging yields an (m,n) system of G* = Al-2 type 

L-2 

m a + n a = ±N6(a,l) + ±J2 A *,bm b , a = 1, 2, . . . , L - 2 (2.3) 

b=i 

where 

AT - 2m - 4n 2 - 6n 3 2(L - 2)n L _ 2 

= 2(1311 ( 2 - 4 ) 

and 

N ^ mi ^ m 2 ^ ••• ^ m L _ 2 ^ (2.5) 

Here A is the A^_2 adjacency matrix. Notice that N and m a must all be even. 

Starting with the (m, n) system, it is possible to eliminate the particle numbers n a in favour 
of the dual-particle numbers m a or vice-versa. It is this fact that underlies the duality between 
particles and dual-particles. Which is to be regarded as the more fundamental is a matter of choice. 
Explicitly, using column vectors, the (m, n) system can be written in the alternative forms 

m + n = \{Nei + Am), n = ±(Ne 1 -Cm), m = C' 1 (N e x - 2n) (2.6) 

where 

m = (mi,m 2 , . . . ,m L _ 2 ) T , n = (n 1 , n 2 , . . . , n L _ 2 ) T , e 1 = (1, 0, 0, . . . , 0) T (2.7) 
and the inverse Cartan matrix of Al-2 is 

C ab = J^i min ( fl ( L - 1 - 6), 6(L - 1 - a)) (2.8) 

Tl'- Let L = \_N/2\ + 1 and let uq be the number of tower particles. Then the geometric 
packing constraints to accommodate n a particles of type a and m a dual-particles of type a along a 
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path of N steps on G = Tl are 

N = n + 2ni + 4n 2 + 6n 3 + 8n 4 H + 2(L - l)n L -\ 

mi = no + 2n 2 + 4n 3 + --- + 2(L-3)n L _2 + 2(L-2)n L _i 

^ (2.9) 
= + 2n L _ 2 + 4n L _i 

m L _ 2 = no + 2n L _i 

mi-i = n 

Rearranging yields an (m, n) system of G* = T' L-1 type 

L-l 

m a + n a = ±iV£(a,l) + ±^A a)6 m&, a = 1, 2, . . . , L - 1 (2.10) 

6=1 

where ^4 is the T' L1 adjacency matrix. Explicitly this yields the inversion of (|2.9p in the form 

n\ = 7}N + \rri2 — m\ 

n a = \m a -\ + \m a+ i - m a , 2<a<L-2 (2.11) 

Notice that for each a, m a = N mod 2 and 

N ^ mi ^ m 2 ^ • • • > m L _i ^ (2.12) 
This (m, n) system can also be written in the alternative forms (|2.6f) where 

m = (mi,m 2 ,... ) m i _i) T , n = (m, n 2 , • • • , n L -i) T , e\ = (1, 0, 0, . . . , 0) T (2.13) 
and, explicitly, 

C-^mi^a,^, C- 1 e 1 = e = (l,l,...,l) T (2.14) 
2.4 Decomposition of paths into particles 

In addition to translational motion, pure particles can be cut into many pieces (with other particles 
inserted) and turned upside-down. To see the particle content of an arbitrary path, as in Figure 3, 
we need an algorithm to decompose the path into pure particles and complexes of overlapping 
particles. Since the particles are composites, the dual-particles must be carried along with the 
particles. Moreover, the algorithm is conjectured to yield a bijection so that it is possible to return 
to the original path given the configuration of particles and dual-particles. 

The particle decomposition algorithm we use is based on the bijective decomposition of War- 
naar [381 139] but is modified to be more symmetric. This modification is required to make contact 
with the various lattice models in the next section. Explicitly, the algorithm to decompose an arbi- 
trary Al or Tl path to reveal the particle content and configurations of particles and dual-particles 
is as follows: 

1. Identify any flat segments corresponding to tower particles and decorate them with dual- 
particles. These tower particles automatically separate the path into complexes with respect 
to the initial baseline at height 1. A complex with respect to the current baseline is any 
segment of the path, starting and ending at the baseline, that is not a pyramid. 



S 



1 

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 



Figure 3: An arbitrary path decomposes into pure particles and complexes of overlapping particles. 
In this way pure particles not only move (have different locations) but also can be cut apart and 
turned upside-down. Even though a particle can be cut into many pieces with other particles 
inserted, the area (shown shaded) of the pure particles (pyramids) is preserved along with the 
decorating dual-particles. In this example there are particles of type 1, 1, 2, 2, 3 at j = 1, 5, 5, 13, 13 
and dual-particles of type 1 at j = 3,7,9,10,12,14,16,17 and type 2 at j = 10,16. Note that 
particles (dual-particles) of different types can occupy the same lattice site. 

2. For each current baseline, separate the pure particles (pyramids with respect to the current 
baseline) from the complexes and decorate these pyramids with particles and dual-particles. 

3. If there are no remaining complexes then terminate the algorithm. For each remaining com- 
plex with respect to the current baseline, identify the left-most and right-most global maxima 
and connect these with a new baseline. The left-most and right-most global maxima may 
coincide and in this case no new baseline is drawn. Decorate the new baseline with a particle 
placed at the midpoint with a type appropriate to its height above the previous baseline. 

4. From each left (right) maxima, follow the profile of the complex moving continuously down 
and to the left (right), and, where a valley is encountered, draw a further baseline at the height 
of the highest peak neighbouring the valley. This process thus introduces new baselines at the 
heights of secondary peaks as they are encountered. Decorate all sloping segments (not valleys 
or peaks) of the complex, outside or at the endpoints of the new baselines, with dual-particles. 

5. For each of the new baselines generated reflect the path constrained within it about that 
baseline and identify each as a current baseline. Go to step 2. 

The iterative stages of this particle decomposition algorithm are illustrated in Figure 4. 

The pure particles are solitons in the sense that they can move together from far apart, interact 
and then re-emerge as unchanged pure particles. This holds generally for an arbitrary number of 
particles but is illustrated for particles of types 1 and 3 and particles of types 2 and 3 moving 
through each other in Figures 5 and 6 respectively. 

2.5 Particle content and string patterns 

The content of particles and dual particles is conveniently encoded in a pattern of 1- and 2-strings 
in a number of physical strips as in Figures 5 and 6. The strips can be displayed horizontally or 
vertically as convenient. Anticipating the connection with lattice models, we call such patterns 
string patterns or zero patterns. A particle of type a at position j corresponds to a 2-string in strip 
a at position j, that is, a pair of zeros at the two edges of the strip with a common coordinate 
along the strip. A dual-particle of type a at position j corresponds to a 1-string in the center of 
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Figure 4: Illustrative example showing three iterative stages of the algorithm for decomposing a 
particular Al path to reveal the particle content and configurations of particles and dual-particles. 
In the first stage, a new baseline is drawn at height 6. The profile is traced down and to the left 
and down and to the right. This introduces (upside down) particles of type 1 at positions 5 and 27 
and the indicated decorations. In the second stage, a new baseline is drawn at height 3. The profile 
of the complex between sites 9 and 23 is traced out adding the indicated dual-particle decorations. 
In the third stage, a new baseline is added at height 5. The profile is traced out and the final 
decorations added. 
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Figure 5: A particle of type 1 moving from right to left through a particle of type 3. The con- 
figurations as the particle of type 1 progresses further to the left are obtained by reflecting in the 
vertical. Like solitons, after the interaction, the particles maintain their integrity. On the right, 
the type 1 and 2 particle and dual-particle content is shown as a string pattern in two strips. 




123456789 10 




123456789 10 




123456789 10 



Figure 6: A particle of type 2 moving through a particle of type 3. Like solitons, after the interac- 
tion, the particles maintain their integrity. On the right, the type 1 and 2 particle and dual-particle 
content is shown string pattern in two strips. 
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strip a at position j. Although the positions j are used to initially generate the string patterns, 
only the relative ordering of 1- and 2-strings in each separate strip is relevant. String patterns with 
the same relative orderings of 1- and 2-strings are regarded as equivalent. 

We assert, without proof, that the particle decomposition provides a bijection between RSOS 
paths and distinct string patterns. The RSOS paths and string patterns are known to give the 
same Unitized characters for arbitrary system size N. In addition, we have checked this conjecture 
explicitly for sizes out to N = 16 by coding the decomposition algorithm in Mathematica [45] . We 
hope to give a proper proof of this conjectured bijection elsewhere. 



2.6 (/-binomials and particle energies 



The distinct string patterns in strip a with n a 2-strings and m a 1-strings are enumerated com- 
binatorially by a (/-binomial generating function [ m "^ na ] q as illustrated in Figure 7. The strip 
configurations are assigned a weight q E and graded by the energy E. The lowest energy config- 
uration, with all 1-strings above all 2-strings, is assigned the relative energy E = 0. Each time 
a 1-string is pushed down through a 2-string, the energy increases by one unit. This excitation 
energy is associated with the 1-strings or dual-particles. 



1 + q + 2q 2 + q 3 + q 4 

Figure 7: Enumeration, by the (/-binomial L] = 1 + q + 2q 2 + + (/ 4 , of distinct string patterns in 
strip a with m a = n a = 2. For convenience, the strips are displayed vertically with position j = 1 
at the top. 



Integer string quantum numbers 
uniquely label the energy levels of dual-particles in each strip. These are defined by 



r(«) 



{number of 2-strings above the 1-string labelled j in strip a} 



(2.15) 



(2.16) 



Clearly, 



(2.17) 
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In addition to the excitation energies described by g-binomials, there is an energy associated 
with the creation of the particle content, that is, particles or equivalently dual-particles. In terms 
of dual-particles, this energy is given by the Cartan matrix. The total energy assigned to a string 
pattern is 

a ja 

where the first term involves the central charge c given by (jl.2p or c = 1, the second term is the 
creation energy of the particles and the third term is the excitation energies arising from the various 
strips as in Figure 7. The central charge term is added to make contact with unitized characters 
given by the generating function or configurational sum 

config 

where the configurations are enumerated by paths or string patterns. The total energies (12.181) 
including the central charge term are correctly recovered in the TBA calculations of Sections 3 and 
4 based on the string patterns. 



2.7 Ai finitized fermionic characters 

In the Al case, the generating function for the spectrum of energy (I2.18|) in the vacuum sector 

is [HE] 

xS(*) = t c/24 E ^ TCm ffp m +n 1 (2 - 20) 

(m,n) a—l * 

where C is the Al-2 Cartan matrix. This is precisely the finitized vacuum fermionic character of 
the unitary minimal model M(L,L + 1) with central charge c given by (jl.2p . In the limit N — > oo, 
with N even, one obtains the vacuum fermionic character 



„t m T Cm L 2 . 1 ( , \ 

c/24 \ ^ TT \2K m a-l + m a+l) 



m a =0 ( 

where tul-i = and the q- factorial is 



.„„-,., 1---, ( Q ) mi a=2 m ° 

in 2 ^ . . . ^ rn L — 2 ^® 



(2.21) 



(q)m = U(l-q j ) (2.22) 

i=i 

2.8 Tl finitized fermionic characters 

Combinatorially, in the Tl case corresponding to the XXX model with central charge c = 1, the 
finitized Tl vacuum character admits several different forms related to spinons, Young tableaux 
and rigged configurations. We recall here the finitized [IB] and infinite N forms [471 14*8] of this 
vacuum character. 
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JV 

N/2 - S z 



In the spinon formulation, the Unitized Tl vacuum character admits the fermionic form 

Setting q = z = 1 gives the counting of states 

X N )> 



N/2 

r l/24 £ , 
S z =-N/2 



N/2 
5 z =-iV/2 



JV 

A/2 - 5 2 



In the XXX spin chain language, JV + , A_ and S z are quantum numbers. A + is the number of up 
spins and A_ is the number of down spins so that the number of spins is JV = N+ + A_ and the 
^-component of spin is S z = k(N+ — A_). Defining the quantum-dimension 



Xn{z) 



z — z 



— , n G Z 



(2.23) 



there is also a bosonic formula for this character [47j 

XiV-2iV + +l(z)| 



LiV/2j 

-i/24 

7V+=0 



/yiV/2-iV + ) 2 


- AT " 




JV+- 



(7V/2-7V++1) 2 



LJVi - 1 



(2.24) 



As a consequence we have 



L7V/2J 



2 Ar = X& ) (M) = E (^-2A+ + 1)Z(A,A + ), Z(A,A 4 

jV + =0 



JV 
JV, 



JV 
JV, - 1 



(2.25) 



and JV — 2JV_|_ + 1 = 25 + 1 is the degeneracy of the spin sector with total spin S. 

An alternative fermionic representation of the Tl vacuum character is given by |47j 



" 1/24 E ^™ W 2j+l( Z ) 



L^v/2J 

,hm T Cm TT f™a + n Q 



(m,n) 



n[' 

0=1 



(2.26) 



where C is the Ti^ Cartan matrix. Let us consider the projection of this character onto the z° 
s£(2) charge sector (S z = 0) 



Xo?(<?) 



-i/24 

(m,n) 

-1/24 £ g 

(m,n) 



74 



m Cm 



^ m T Cm 



[_N/2\ 

n 

o=l 

n[ 

a=l 



rm a + n a 




j-l/24 


- JV - 


V m a . 


g 




.A/2. 



(2.27) 
(2.28) 



This is convenient since, from the lattice, this projected character contains all of the distinct 
eigenvalues (as Laurent polynomials) and removes the degeneracies. In the second formula we have 
taken an arbitrary integer K > [-V/2J = L — 1 and extended the definitions of m and n by 



n a = 0, m a = m L -i, a>L-l 



(2.29) 
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with C = T' K . Taking the limit K — > oo followed by N — > oo, with N even, we then obtain the 
projected vacuum fermionic character 

r,h rn T Cm . 1 I ^ _i_ ™ \ -, 

m a — (mod 2) v ' - 1 a=2 y 
^ mo ^3 . . . 

where C is now the Cartan matrix of = and the sum is over asymptotically constant 
configurations m for which there exists some finite K > such that m a = tuk whenever a > K. 

2.9 One-dimensional configurational sums 

Typically, unitized characters can be expressed as one-dimensional configurational sums, that is, 
via an energy statistic applied directly to RSOS paths. These one-dimensional configurational sums 
first appeared in the work of Baxter [3] but in the context of off-critical Corner Transfer Matrices 
(CTMs) where the elliptic nome q (|1.3p is related to the departure from criticality. By contrast, 
we work here only at criticality and q is always the modular nome. 
The one-dimensional configurational sums are defined by 

X^ ) c {q) = Y,Q lT '^ ima3 - 1 ' (T ^ a3+1 \ *o = a, a N = b, a N+1 = c (2.31) 
M 

where H{aj-\, crj, is a local energy (energy density) and the sum is over all RSOS paths 

a = {o"o, <Ti, . . . , Ojv+i} on Al or Tl with Oj €{1,2,..., L}. These one-dimensional configurational 
sums satisfy mathematically powerful recursion relations in N relating different boundary conditions 
specified by a, b and c. In the vacuum sector with N even, a = b = 1 and c = 2. 

The form of the local energy function is not uniquely determined due to the possibility of local 
gauge transformations. Let us choose, however, the local energy function to be of the simple form 



1, (aj-i,crj,aj + i) = (1,1, 1) 

i, (aj-i,aj,a- j+1 ) = (1,1,2) or (2,1,1) 

1, \<Tj-l - <7j + l| = 2 

0, otherwise 



(2.32) 



It is then confirmed that, for both Al and Tl paths, the energy statistic 

1 N 

E = ~^ + 2^ j H ^-i> ( 2 - 33 ) 

3=1 

exactly reproduces the string energies (|2.18p and hence 

X { %) = q - c/U xS 2 \q) (2.34) 

This relation has been checked directly in Mathematica [H] for all sizes up to to N = 26 and 
confirms the equivalence of (|2.18p and (I2.33p . This largest size involves 10.6 million paths. 
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2.10 Combinatorial classification of states 



The combinatorial classification of states for the A§ minimal model for paths of length N = 8 is 
shown in Figures 8 and 9. In this case there are 14 such states with the unitized vacuum character 



The combinatorial classification of states for the T5 XXX model for paths of length N = 6 is shown 
in Figures 10 through 12. In this case there are 20 such states with the unitized vacuum character 



In both of these cases, there are 3 physical strips and we see that the patterns of 1-strings and 
2-strings in these strips arising from the particle decompositions of the paths exactly match the 
patterns of eigenvalue zeros in the upper-half u-plane obtained by direct diagonalization of the 
commuting double row transfer matrices with N faces along a row. 

The particular system sizes in these examples were chosen for convenience in illustrating the 
combinatorial classification of states. Many values for L and N have been checked, in particular 
considerably larger system sizes, and in all cases perfect agreement was found between the combi- 
natorial picture in terms of paths and the numerical location of zeros of eigenvalues of the transfer 
matrices. 

2.11 Bethe roots and transfer matrix 2-strings for XXX model 

In Figures 10 through 12, relevant to the XXX model, we have also shown the eigenvalue zeros of 
Baxter's auxiliary matrix Q(u) as in (|4.3|) and (|4.4|) . These zeros are precisely the Bethe ansatz 
roots. The transfer matrices D q {u) are completely determined by knowledge of the matrix Q(u) 
through Baxter's T-Q relation (|4.3j> and the fusion hierarchy. 

It is readily confirmed, in this example, that an eigenvalue Q(u) has a j-string at the same 
relative position as each 2-string in strip j of the transfer matrix eigenvalues with j = 1,2,... , N/2. 
Identifying S = no/2 as the total spin, we see that there are 5,9,5, 1 eigenvalues with total spin 
S = 0,1,2,3 respectively in the S z = sector. The zero content of an eigenvalue Q(u) thus 
determines the complete particle content n a with a = 0, 1, ... , N/2 — 1. The dual-particle content 
then follows from the (m,n) system (|2.10p . Again we have verified these observations for many 
different system sizes N. We further conjecture that these observations hold true for all N and 
that in the limit N — > 00 the locations of the j'-strings of Q(u) exactly coincide with the locations 
of the 2-strings of the transfer matrices in strip j. 



The bulk face and boundary triangle Boltzmann weights of the critical unitary minimal RSOS 
models [251 ES] on Al in the vacuum sector are 



X ( % 5 (q) = q- 1/30 (l +q 2 + q 3 + 2q A + q 5 + 2q 6 + q 7 + 2q 8 + q 9 + q 10 + q 12 ) 



(2.35) 



X ( %(q) = 9~ 1/24 (1 + q + 2q 2 + 3q 3 + 3g 4 + 3q 5 + 3q 6 + 2q 7 + q 8 + q 9 ) 



(2.36) 



3 Unitary Minimal RSOS Lattice Models 




(3.1) 
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Figure 8: Eigenvalues 1 to 7 for A§ with N = 8 showing the energies (omitting the — c/24 term), 
paths, predicted patterns of zeros and actual patterns of zeros (rotated through 90 degrees) from 
numerical diagonalization of transfer matrices. Only the relative order of 1- and 2-strings in each 
strip is relevant. The relative order across strips is not relevant. 
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Figure 9: This is a continuation of Figure 8 showing eigenvalues 8 to 14 for A§ with N = 8. 
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Figure 10: Eigenvalues 1 to 7 for the XXX model with N = 6 showing the energies (omitting 
the — c/24 term), paths and actual patterns of zeros (rotated through 90 degrees) from numerical 
diagonalization of transfer matrices. Only the relative order of 1- and 2-strings in each strip is 
relevant. The last column shows the numerical patterns of zeros for the eigenvalues of the auxiliary 
matrix Q. The Q eigenvalue has a -/-string at the same relative position for each 2-string in strip 
j of the transfer matrix eigenvalues. The path out to j = 7 is needed to apply (2.33). 
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Figure 11: This is a continuation of Figure 10 showing eigenvalues 8 to 14 for the XXX model with 
N = 6. 
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Figure 12: This is a continuation of Figures 10 and 11 showing eigenvalues 15 to 20 for the XXX 
model with N = 6. 
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where 



s(u) = — - T = r , ^ = e» x = e* A , A= — — , S a = s(a\) (3.2) 

sin Ax — x 1 L + 1 

Commuting double row transfer matrices can now be defined following |35j . 
3.1 Bethe ansatz 

The double row transfer matrices of the unitary minimal models in the vacuum sector satisfy the 
T-Q functional equation [10] modified to account for the boundary 

s{2u) D{u + ^)Q(u) = s{2u + X)f(u + ^)Q(u - A) + s(2u - X)f(u - -)Q(u + A) (3.3) 

where f(u) = s(u) 2N and Q(u) is an auxiliary family of matrices satisfying [Q(u),Q(v)] = 
[Q(u), D(v)] = 0. The same functional equation is satisfied by the eigenvalues D(u) and Q(u) 
and so the Bethe ansatz equations result by setting 

n 

Q( u ) = s ( n — Uj)s{u + Uj) (3.4) 

3=1 



where Uj are the r = 2n Bethe roots. The physical analyticity strip of D(u) is —A/2 < Ren < 3A/2 
the physical analytici 
Let D(u) = Dl and 



3 

and the physical analyticity strip of D(u) = D(u + |) is — A < Ren < A. 



D\ = D q (u + kX), Q k = Q(u + kX), s k (u) = s(2u + kX), f k (u) = (-l) N s(u + kX) 2N (3.5) 

then the T-Q relation implies that the eigenvalues D{u) are determined by the eigenvalues Q{u) in 
the compact form 

~ si/ 1/2 Q-i + s_i/_ 1/2 Qi 

3.2 Functional equations 

The fused transfer matrices D q satisfy the fusion hierarchy of functional equations [30L [35] 

S q - 2 S2q-lD q D 1 q = S q - 3 S 2 qf q D q ~ 1 + Sg_lS 2 g-2/g-l£ , o +1 , q = 1,2, . . . , L — 1 (3.7) 

subject to the closure conditions 

£>o 1 = 0> ^0 = /-!^ ^>o=0 (3.8) 

Starting with the fusion hierarchy induction can be used to derive the T-system of functional 
equations [32j [35] 



s q - 2 s q D g Dl = s. 2 s 2q f^ l f q I + s 2 q _ l Dl +1 D\'\ q=l,2,...,L-l (3.9) 
For q = 1, this is just the usual inversion identity given by the fusion hierarchy with q = 1. 
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If we further define 

d q = q ~ l \ ? , 9 = 1,2,. ..,£-2 (3.10) 
then the inversion identity hierarchy can be recast in the form of the Y-system 



S-2SoS2qS2q+2f-lfofqfq+l 

^t 2 D\ Yr , slD\Dr\ . 



i + ; u i + ; 1 / j = (i + dr i + dp) (3.ii) 



with closure 



dg = d^" 1 = (3.12) 

We wish to solve the Y-system (13. lip subject to the constraints of periodicity, conjugation and 
crossing symmetries, analyticity and asymptotic limits. The periodicity, conjugation and crossing 
symmetries are 

D q {u) = D q (u + w), D q (u) = D q {\-u) (3.13) 



Dl(u) =D q {u) = D q ((2-q)X-u), d«(u) = d q (u) = d q ({l - q)X - u) (3.14) 

To discuss analyticity, we consider the eigenvalues of the transfer matrices D q (u) and d q (u) at each 
fusion level in their respective analyticity strips 

-|A<Reu<-^A, -^-X<Reu< A (3.15) 

Defining the shifted transfer matrices 

D q {u) = D q (u + ^^-X), d q (u) = d q (u + ^-^-X) (3.16) 



it follows that these transfer matrices have the common analyticity strip 

- A < Re(u) < A (3.17) 
and satisfy the same periodicity (|3.13p with crossing symmetries 

D q {u) = D q {-u), d q (u) = d q {-u) (3.18) 
In terms of shifted transfer matrices, the TBA functional equations take the form 

d q {u - ^) d q { u + ^) = (i + d q -\ u )) (i + d q+ \u)) (3.i9) 

Lastly, the asymptotic values d 3 (+ioo) were computed in [32] 

sin[j<9] sin[(j + 2)9] sin 2 (j + 1)0 1 a _ stt 
sin 2 9 sin 2 

Here, s is a Kac label and plays the role of selecting the eigenstates of the transfer matrix that, 
in the scaling limit, will be in the (r, s) sector of the conformal field theory. More precisely, s is a 
good quantum number for the finite system corresponding to the braid limit u — » ioo. In contrast, 
r is only a good quantum number in the limit N — > oo and so is not accessible directly on a finite 
lattice (see after (|3.94p ). In the vacuum sector of interest here, (r, s) = (1, 1). 



d3 (+^) = -^-dri — — = l* a e = TTT (3 ' 20) 
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3.3 Analyticity and solution of TBA 

We use analyticity properties of the transfer matrices to transform the Y-system ()3.11j) into the 
form of Thermodynamic Bethe Ansatz (TBA) equations. For the Al lattice model, these TBA 
equations take the form of coupled integral equations whose structure is governed by the Dynkin 
diagram G* = A^_2- The analyticity is encoded in the relative locations of the 1- and 2-strings in 
the analyticity strips (|3.15j) . The TBA equations can be completely solved in the scaling limit for 
the locations of the 1-strings as well as the excitation energies. Following [37], the strategy is as 
follows: 



1. The special case q = 1 of (13. 9p provides the expression for the energies or, equivalently, the 
eigenvalues of the transfer matrix. This expression contains an unknown function (pseudoen- 
ergy in the first strip) and a number of unknown real numbers (location of 1-strings in the 
first strip). 

2. The Y-system (|3.1ip can be transformed into the set of TBA equations. 

3. The Y-system (|3.1ip allows us to access the location of the 1-strings by a set of coupled 
auxiliary equations, that perform the job of quantisation conditions. 

4. The simultaneous solution of TBA plus auxiliary equations completely fixes the pseudoener- 
gies, the locations of the 1-strings and the energy expression. 

3.3.1 Analyticity: the energy 

We aim at finding the transfer matrix eigenvalues in the physical strip ^ < Reu < 4p. We start 
from the case q = 1 of (|3.9|) , written for the eigenvalues of the corresponding matrices 

/sin2A\ 2 s-xai ^ D l (u + A) = 1 + d\u) (3.21) 
V sin A J s_ 2 s 2 

Ignoring the term d}{u) which is exponentially small for large N, this functional equation factorizes 
into bulk and boundary equations with D l {u) = ftbuik( w ) K o( M ) 



[«bulk(«) «buik(« + A)] = f-if_ 



sin(A — u) sin(A + u) 1 " ' 



sin 2 A 



Reu€(-^) (3.22) 



i , , . S-oSo ( sin A \ 2 sin(A — 2u) sin(A + 2u) ( sinA\ 2 „ , A A. 

K (u )kq(u + X = — — =-7^7 . — ^TT , RfiWe (--,-) 

s_isi \sm2Ay sm(2A — 2u) sm(2A + 2u) \sm2Ay 2 2 

(3.23) 

The expressions for these factors are given in [19] and immediately extend to the whole physical 



strip Reu G (— -|, |A). It follows that the "finite" part of the transfer matrix is defined 

""'"^ [JySI^ ^"^ 1 ' B"«(4§A). (3.24) 



Here we do not make use of these factors we only need to know that, from [49], they do not introduce 
new zeros or poles to Di. t Ju) in the physical strip. It follows that the finite transfer matrix has 
precisely the same zeros as D l (u) and no poles. The functional equation (|3.2ip now takes the form 

^finit>-^) bl^ e {u + ^) = l + d\u), ReuG(-^). (3.25) 
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Let y^ denote the location of the zeros of L' 1 (tt) such that 

D^(±%P) = D\±%^) = 0, yW > 0. (3.26) 



7T 7T 

1/ 



To remove these zeros from D (u), we construct the function 



mj 



z(«) = n 



fc=i 



Re« G (—A, A) (3.27) 



In the indicated region this function has the same zeros as D (u) and no other zeros or poles. 
Moreover, it satisfies the inversion identity 

Z(u-|)Z(u + = (3.28) 



Let us define the function 



A (x) = gfefc , . R(.»G(-A.A) (3.29) 



Z(u) 



■£+1 



which is free of zeros and poles in the analyticity strip Imx G (— ir, 7r). Dividing (|3.25p by (|3.28j) 
we then obtain 

A(x + i l) A (x-^) = l + e\x) , Imx 6 (-- |) (3.30) 



where 



&(x) = d?(u)\_ ta . (3.31) 

IU— L+l 



The left side of (|3.30p is free of zeros and poles in the indicated region so the right side never 
vanishes. The function d l (u) vanishes at u = by (|3.10p and (|3.8p . We conclude that 

l + e 1 (x)>0 for x€K. (3.32) 

so that we can take the logarithmic derivative 

— log^l(x + y) + — log^l(x-y) = — log(l + e 1 ( 2 ;)) (3.33) 
Taking Fourier transforms 

Hf,k) = — dxf(x)e- ikx , f(x)= dkF(f,k)e ikx (3.34) 

^ J — oo J — oo 

we obtain 

F(—logA(x),k) 2cosh(fc~) =^f-^-log(l + e 1 (x)),A; N ) (3.35) 
Vox / 2 Vax / 

Transforming back, we obtain a convolution with a kernel and an integration constant 

7, , + const = (K * log(l + e 1 )) (x) + const (3.36) 
2vrcosh(x-y) 



where 



(/ * = ^ [°° dy f{x - y) g(y), K(x) = —L- . (3.37) 
27r /_„ cosnx 
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Using (|3.24p , (|3.29p and the needed functions we can reconstruct the initial double row transfer 
matrix eigenvalues log -+- ^ct)- Actually, here we are concerned with the finite size corrections 
to it namely those terms that characterise the scaling limit of the theory. These terms behave as 
1/N in the log of the transfer matrix. The bulk term in (|3.24p is of order N; the boundary one and 
the integration constant are of order 1 and will be ignored. What remains expresses the scaling 
part of the free energy (in (|3.36[) we ignored the constant) 



_,i ,A ix , 



/seal. (•'') ~ 

(K*log(l+e 1 ))(x)-^ 



logA{x) - Z( 



L + V 



(3.38) 



mi 



k=l 



(1) (1) , 

vi, — x yl + x 
log tanh — h log tanh — 



We can fix x = namely choose the isotropic point of the system and get 

mi 



/seal. = - log £>fi 



A , 



finite V 



2tt 



" d ,!2^±£M_ 2E l„ g taoh 

k=l 

1 



2ir cosh(y) 



(i) 



(3.39) 



— E + higher order corrections in — 
N N 



and E is the energy in the language of conformal field theory, corresponding to (|2. 18j) 

c 



E 
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+ conformal dimension 



(3.40) 



So far, the zeros' positions and the function e 1 are unknown: we will evaluate them in the next 
sections. 

To summarise, starting from (|3.21|) . we have discussed the technique, which we will largely 
use later, consisting in the removal of all zeros and poles from the analyticity strip, then using 
the Fourier transforms on the non-singular part of the functional equation to extract the unknown 
function £) 1 (it). 

3.3.2 Analyticity: the TBA 



We aim at solving the Y-system (I3.19p . for q generic. We start examining zeros and poles of d q (u) 
in the analyticity strip (I3.17p . using the definition (|3.10p converted to the tilde form 



<P(u) 



sin 2 (2ii) 



D q ~ 1 {u)D q+1 { 



u 



sin(2n - (q + 1)A) sin(2n + (g + 1)A) 

The bulk term in square brackets vanishes if 

u = ± — - — A 



sin 2 A 



-i 2N 



(3.41) 



(3.42) 



These points are outside or on the border of the analyticity strip for all q so they will play no role 
later. On the contrary, for the special case q = 1, the transfer matrix of level zero has a zero at 
u = 



D°(u) 



-1 



sinti 
sin A 



2N 
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that implies the same type of zero for d}{u) so we need to remove it from the case q = 1 in (|3.19j) : 
to that purpose we need a function with a zero of order 2N in u = 0, with no other poles or zeros 
in the analyticity strip and satisfying 



i 



We notice that iV is even, here. Clearly we have (an overall sign is irrelevant) 

'L + l 



f(u) = i tan 



-u 



(3.43) 



(3.44) 



for q = 1 only. The first factor in (|3.4ip has poles at the points indicated in (|3.42p . that lie outside 
the analyticity strip; it has also a double zero in u = that we will remove as in (|3.43p . (I3.44p with 
the same functions 



g{u) = [f(u)f 



lan 2 [ £±iu 



9{u 



\)g(u+±) = l 



(3.45) 
(3.46) 



In addition, we have to remove the zeros of the strips q ± 1. These zeros are in the center of 
the analyticity strip so, remembering (|3.27p and (|3.28|) . we introduce 



L-2 rrij 

^) = n n 

j=i k=i L 



tanh-^ + ^-^ tanh-^ + 1) ^ + ^ 



Q) 



(3.47) 



(3.48) 



where the adjacency matrix (12. ip selects the appropriate strips. As indicated, this function satisfies 
an inversion identity. Now we divide the functional equation (13.19)) by (|3.43p . (I3.46p . (13.48P and 
obtain a newer equation for the analytic and non-zero function A q (x) 



A q (x) 



A q (x + ^)A q (x-^) 



d q (u) 



[f(u)}^^g(u)Zi(u) 



Re u G (-A, A) 



■ L+l 



(l + £ q - 1 (x))(l+i q+l (x)), Imxe(- 



7T 7I\ 

2' 2' 



logd^fa)) ( „ 
2tt cosn(x — y) 



As in (|3.30l and I3.36H . we solve for A q (x) and obtain 

L-2 

logA q (x) = J2A qJ 
i=i 

L-2 

= ^A 9j (if*log(l + ^))(z) + C^. 
We now replace the definition of A q (x) and get the TBA equations 



(3.49) 
(3.50) 

(3.51) 



loge«(x) = 5i,,JVlog[/(- 



IX 



'L + l 

L-2 rrij 

+ ^g n n [ tanh 

j=l fc=l 



+ \ogg 



L + l 



+ 



(3.52) 



(i) 

1 »>' .an!,:'' • ^ 1 ' 



L-2 



+ J2A q , j {K*log(l + ^))(x) 
i=i 



1 With more general boundary conditions, these two functions need not to coincide and g must contain the boundary 
parameter £. 
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The +ioo asymptotic values (|3.20p fix the integration constants; we make use of the integral 

dy — = 7r (3.53) 

3 cosh y 

and also observe that (|3.49|) being a non-zero function, the expression 1 + cannot vanish on the 
real axis and its sign can be evaluated from its asymptotic value (|3.20p 

sin 2 (o + 1)A 

l + e q (+oo) = ; >0 (3.54) 

sin A 

therefore on the whole real axis we have 

l + e q (x)>0 i£K, q = l,...,L-2 (3.55) 

and its logarithm does not introduce imaginary terms. 
We can now evaluate the limit in (13.521) 



loge*(+oo) = \Ys A i,i Ml + ^(+oo)) + (3.56) 

where logarithms are taken in the fundamental branch. Using (|3.20p . (|3.55p we conclude that 

C (9) =0 q = l,...,L-2. (3.57) 

Here we restrict ourselves to the vacuum sector but in other cases the integration constants could 
be non-zero, usually multiples of i tt. The L — 2 TBA equations are now 

log e q (x) = S 1>q N log tanh 2 | + log tanh 2 | + (3.58) 



L-2 rrij (j) (j) L-2 

log I] II [tanh X -^- tanh ' J + £ A g , (K * log(l + e*)) ( 



3.3.3 Analyticity: the zeros 

The zero's positions are still undetermined and, again, the functional equations (|3.19p will help us. 
We shift them of ±^ so we can have one of the two forms 

$4(u)$(u + \) = (l + d q - 1 (u+^))(l + d q+1 (u+^)) (3.59) 



2 / V v 2' 

d q {u-\)d q {u) = (l + d q - l (u-^(l + d q+1 (u-^ (3.60) 

and if we are on a zero of d q {u) both expressions vanish. The zeros of d q (u) are in both the strips 

q ± 1 (|3.4ip . Let's take a zero of q + 1, no = ' ; in that case also d Q+2 (uo) will vanish but not 
d q+A {u ) so 



d"(u o ) = => (l + ^- 1 (n + ^))(l + ^ +1 (n + ^)) =0 
d q+2 (u ) = => ^ + d q+1 (u + ^)){l + d q+3 {u + ^))=0 (3.61) 
d q+ \u )^0 =► ^ + d^(u + ^){l + d q+5 (u + ^))^0. 
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The non- vanishing of 1 + d q+3 (uo + ^) forces 1 + d q+1 (uo + ^) to vanish on the zeros of strip q + 1 
and to be different from zero on the zeros of the remaining strips. The same is true for the 'minus' 
form in (I3.60p so we state that 

= 1 + * ± |) = 1 + e q T q) (3.62) 

Actually, by (|3.41|) . d q (u) also vanishes if u is a 2-string of q + 1 (and also of q — 1). The 2-strings 

are located in uq = ±A + tw ^ +1 with real w^ +1 \ Taking the "— " case, we use (|3.60p and repeat 
precisely the same construction as in (|3.61|) ; this yields (q + 1 — > q) 

= 1 + ^(^ + ^)=l + ^(-A + ^ + ^)=l + ^(4 9) +4) (3-63) 

This calculation holds true for the other sign so we conclude that a unique equation, equivalently 
in the 'minus' or the 'plus' form, 



= 1 + £*(«=f4) , ve{y^\k = l,...,m q ; w ( h g \h = l,...,n q } (3.64) 



.IT 

2 

fixes both the one- and the 2-strings of the given q strip. We don't need the 2-strings positions but 
we need their order with respect to the 1-strings. Focusing on the "minus" form, we introduce the 
function 

^ q (x) = iloge q (x - i-J = i<5i >g iV log tanh 2 (- - i-J + ilogtanh 2 - i-J + (3.65) 



L-2 rrij ^ (j) ^ ^ (j) ^ 



+ i ^2 X> [ ^ og tan ^ (~ — ^ *t) + tann 



i=i fc=i 
L ~ 2 log(l+#(y)) 



4 



+ 



/oo 
-oo 



2-7T sinh(x — y) 



that is real on the real axis and becomes an odd multiple of ir on one- and 2-string zeros. Differently 
from (|3.58p . here the sum on different zeros is done after taking the logarithm so we can get an 
extremely useful information by keeping track of the winding. On the zeros' positions the principal 
value is not required as both the numerator and the denominator vanish while it is required for 
generic values of x. We introduce the 1-string quantum numbers by 

y q (y[ q) ) =Trn[ q) n[ q) = 1 mod 2 . (3.66) 

The yu' are always single zeros and also each quantum number uniquely fixes one zero. These 
features have been largely observed in our numerical analysis and also in previous cases as in [50J; 
this is also typical of models solved by Bethe Ansatz. For these reasons, these quantum numbers 
were called "non-degenerate" in [50]. In the framework of ABF models, exceptions are known [50] 
and are related to bulk or boundary renormalisation flows and do not play any role here. Using 
(I3.65P and also (I3.20p we can derive the asymptotic values of the function ^ q . At +oo 1 + d? is 
non-zero for all j so we can safely take the limit of the integral part and get zero 

$9(+oo) = q = l,...,L-2 (3.67) 
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From (|3.20p we find that 1 + dP{— oo) is non-zero for all j > 1 but is exactly zero for j = 1 so the 
evaluation of the integral in (|3.65l) requires special care. Actually, the case j = 1 occurs just once 
for q = 2 so in all the other cases we have 

L-2 

qi(-oo) = 27r{8 hq N + l + J2 A qJ m j) > °> ( 3 - 68 ) 

i=i 

We believe that the divergence of log(l + d 3 {— oo)) does not affect the result so we write 

L-2 

^ 2 (_oo) = 2vr(l + A 2,j m j) > ( 3 - 69 ) 
i=i 

These asymptotic behaviours say that the function is globally decreasing. Moreover, the real 
function 

ilogtanh(| - ijj (3.70) 

is monotonically decreasing. We cannot prove the behaviour of the integral part but it is usually 
subdominant so we conclude that the function fy q (x) is monotonically decreasing. This has been 
widely confirmed numerically; again, minor exceptions are known [SD] and are related to bulk or 
boundary renormalisation flows and do not play any role here. 

The +oo asymptotic fixes the smallest quantum number in (j3.66p : indeed, the largest zero is 

(q) 

Vm q SO 

nnW =^(yW)>^(+oo)=0 (3.71) 



so 

b k+l 



n ( f ) >...>n^>n^ +l >...>n^>l (3.72) 



Given a zero with quantum number ni 9 ' , if a lower 2-string moves up and exchanges its position 
with the 1-string, the position y[ 9 gets smaller and the decreasing monotonicity of ^^(x) implies 
that the quantum number has to grow of two. It cannot grow of one because it must be odd. At 
the same time, the so called non-negative quantum number li 9 ' grows of one by definition (|2.16p 
so the two quantum numbers must be proportional 

njf) = 2I^ ] + integer (3.73) 

Imposing (|3.71l and I3.72p we get 

nf = 1 + 2(4 9) + m q - k) (3.74) 

Both +oo and — oo asymptotics are even multiples of 7r; as one- and 2-strings are quantised 
with odd multiples of it, the number of odd positions from oo) to l I / ' J (+oo) must contain the 
number of one- and 2-strings 2(m q + n q ), the 2 comings from the upper and lower half planes, so 

\rr9('_ 00 N ) _ \T/<?('_|_ 00 N ) 

— K - ; [ ^ ' = S 1>q N + 1 + A i,3 m i > 2 K + « 9 ) (3-75) 

n 3=1 

The 1 is the spurious zero u = of (I3.41|) . 

The next step will be to evaluate the scaling limit of these equations. 
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3.3.4 The scaling limit 

As we are interested in the conformal energy (|3.39p . we need to evaluate 

17(1 



E 



lim < — 



00 dy lo^l+e 1 ^)) 
2vrcosh(y) 



2 log tanh 



k=l 



(3.76) 



so we are pushing the number of sites to infinity. In doing this, we need to pay attention to the 
movement of the one- and 2-strings in consequence of the growth of N. They depend by iV with 
the law 

(3.77) 



yk\ N ) = Vk ) + iogA r >o 



where the scaled position y^ is n °t constrained to be positive, as it was for its lattice partner. We 
also need to consider that the reduced transfer matrix eigenvalues e l {y) also depend by ./V both by 
the TBA equations (|3.58p and by the zero positions. 

The "zeros term" in (I3.76P can be easily worked out using the scaling (I3.77|) 



lim 



N y[ 1} 
•-log tanh 



log N 



1 _ fl W 
—e Vk 

7T 



(3.78) 



The integral term has to be split and treated by changing the integration variables 

log(l + £*(!,)) 



lim 



I 

N 
2^2 



lim 



dy 



00 dy log(l + i^y)) | 







2 cosh(y) 



2tt cosh(y) 

d \og(l +i\y)) 
—00 



(3.79) 



2 cosh(y) 



lim 

N^oo 



lim 

iV^oo 



N 
2^2 

1 

'2^2 



30 d , log(l+e l (y' + logN)) + 
log N 2 cosh (y' + log AT) 

30 . logq+gV + logAQ) 
logiv V eV+N-H-V 



lo * N ff log(l+eV-log^)) 



dy" 



logiV 



2 cosh (y" - logN) 



„ \og(l + e\y" -\ogN)) 

aV N -2 e y" + e -y" 



In the limit, both the integrations go over the whole real axis; the subleading N 2 term disappears 
and we are lead to define the following scaled eigenvalues 



e q (x) = lim s q (x + logN) 

N^oo 

(x) = lim e q {— x — log N) 



(3.80) 



' lower 



The first one comes from the original upper half plane the second one from the lower one. For finite 
N, here we have perfect mirror symmetry between the two half-planes so we don't need to consider 
the lower one. Because of this mirror symmetry, it is convenient to put —x in the definition of the 
lower component. We will show later that scaling the eigenvalues as in (I3.80p is perfectly consistent 
with the TBA and auxiliary equations. 

We finalise the energy expression {y" — » —y") using the mirror symmetry 



1 

'2^2 



/CO f'OO 
dye-y\og{l+e\y))+ / dy e~* log(l + e 1 ^)) 
-00 J —00 

T / dye-y\og(l+e l (y)) 
r J — 00 



(3.81) 
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and finally the conformal energy is given by 

E = -- / dye-y\og(l + e\y)) + 2Y^-e~ti ) (3.82) 

The scaling eigenvalues in (|3.80p can be determined by performing the same limit on the TBA 
equations (|3.58j) . We initially treat the special term occurring for q = 1 

lim dig N log tanh 2 x + l °& N = -^ x . q e~ x (3.83) 

N^oo ' 2 

Then, the limit on the boundary term (second term on the right hand side of (13.580 ) trivially 
vanishes. Among the terms containing the zeros, we look at those where there is a difference 

lim log tanh x + lo ^ N ~ Vk_ = j tanh (3 84) 

while those where there is a sum disappear 

x + log N + y[ q) x + y[ 9) + 2 log N , 
lim log tanh ^- = lim log tanh ^ — = 3.85 

Clearly, these terms were produced by zeros in the lower half plane so their contribution vanishes 
in the upper half plane. The integral term is evaluated using the properties of the convolution 

(K *log(l +e j )){x + logN) = (K * log(l + e j )) (x) (3.86) 

We are ready to recompose the TBA equations 

L-2 ™i X _A3) A . 

loge«(x) = -4 S 1>q e~ x + log ]J ]J [tanh ~^^\ " 

3=1 k=l 

L-2 

+ ^A w -( J pr*log(l + ^))(x) J g = l,...,L-2 (3.87) 

i=i 

Observe that the positivity condition 

l + e q (x)>0 (3.88) 

remains true and the asymptotic value (|3.2U|) is not changed by the scaling limit (|3.8U|) ; that 
asymptotic value was used to fix the integration constants (|3.57j) . We now use the TBA equations 
to evaluate the opposite asymptotic e^ q '(— oo). First, we observe that the zeros contribution is 

lim log T] TT [tanh *"^ 1 " = (3.89) 

j=l fc=l 

because the zeros are in even number, from section [2~3l If we exponentiate (|3.87p . this contributes 
a factor one. The integral term is evaluated as in (|3.53p and following equations, the positivity 
condition (|3.88p ensuring that no imaginary parts are introduced so we have 

lim e q {x)= lim exp(-4 8 X q e~ x ) TT \l + e j (-oo) 2 ^ (3.90) 

3=1 



32 



If q = 1, the first factor forces a transcendentally fast vanishing 

e 9(-oo) = (3.91) 

that is an indirect manifestation of the presence of a sea of 2-strings in the first strip, namely 
pyramids of hight one. For the remaining values of q we write 



L-2 _ j 



Yl [l + e J '(-oo) 2 q ' 3 q = 2,...,L-2 (3.92) 
i=2 

where the range of the product starts from j = 2 because of (|3.91|) . This equation can be solved 
with the technique of [32J: simple trigonometry shows that 

1 + ^-00) = Sin2(g 2 +Q)T (3.93) 
sin t 

is the solution of (|3.92|) for generic q, with two arbitrary parameters a, r to be fixed ( (|3.92p is like 
a second order finite difference equation in the integer q, if the logarithm is taken). The case q = 2 
fixes a = and the case q = L — 2 fixes r = tt/L so we conclude with 

e *(-oo) = ^^-l, q = l,...,L-2, t = j (3.94) 
sin t L 

Actually, the q = 1 case is also contained here. For generic parities of m q , the value of r would have 

been r = rir/L with an integer parameter r in the role of the Kac index of (r, s). The parameter 

s was introduced on the finite lattice in (|3.2U|) ; the parameter r first appears here after the scaling 

limit; in the present case (r, s) = (1, 1) namely we investigate the vacuum sector. 

In literature, the following pseudoenergies are often used 

e q (x) = -loge q {x), q = l,...,L-2. (3.95) 

A standard notation is also the function 

L q (x) = log(l + e 9 (x)) =log(l + e" e<?(:c) ), L,(x)eR for x G R (3.96) 

From (|3.41|) . this function vanishes on the zeros of the adjacent strips 

L q (y < j* ±1) ) = V k. (3.97) 

The quantisation condition (|3.66p . according to the shift of the zeros (I3.77p . becomes 

lim y q (y { k q±1) ) = ^ q (y { k q) + log N) = tt n[ q) n[ q) = 1 mod 2 (3.98) 

so we are lead to scale the function as we did in (|3.8Up . Actually, we recycle all our previous steps 
(I3~m TO . SSm and get 

4> q (x) = lim W(x + logN) (3.99) 

N^oo 

L—2 m r *M _ 

= 45 hq e~ x + iY^ Kr lQ g tanh ( X ~ ) 

r=l k=l 



L -' - log(l + ^(y)) 



dy 



| „_ v 2vrsinh(x-y) 
= 7rn {q) (3.100) 
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Observe that ^f 1 keeps track of the infinite number of 2-strings in the first strip that form the Dirac 
sea. In a familiar terminology, it is the coloured node of the ^4l-2 TBA diagram. Equations (I3.87|) 
and (I3.99P are fully consistent with the scaling limit taken in the energy expression (I3.82|) . 

We can use asymptotics to count the number of one- and 2-strings. Because of (|3.99p . the 
spurious zero u = of (|3.75p now has left the analyticity strips moving toward — oo therefore all 
the odd multiples of tt from if? q (— oo) to ^(-l-oo) are exclusively occupied by one or two strings; 
their amount is m q + n q as the lower half plane is not visible here. The number of odd positions 
from 1 $> q (— oo) to ^ q (+oo) is evaluated from (|3.99p and must be equated to the number of objects 
that can be allocated in the strip, namely m q + n q 

* 9 (-oc) - #«(+oo) l£=? A , 

2^. — - = 222 A ^ m j = m i + n i q = 2,...,L-2. (3.101) 

The case q = 1 simply gives oo = m\ + n\, consistently with the notion that in the first strip there 
is a Dirac sea of n-family particles of high 1, n\ = oo, and mi is not upper bounded. In summary, 
we feed in the functional equations (13.190 the 1-string positions and from the equivalent TBA and 
auxiliary equations we get out the part of the (m, n)-system that survives after taking the scaling 
limit in (1231). 



3.3.5 Energy from excited TBA 

Extending the calculations introduced in [37], the expression for the energy (|3.82p can be cast into 
a closed form. We rewrite the quantization condition in the first strip, as given in (|3.74p with q = 1, 
in the following form: 



2 -e-^ = ^-U\y^) + 2 -e-^ ,3.102) 



We now sum on all the first strip zeros to obtain 

„ mi mi . mi m-2 „(1) „(2) 

f£^" = ^"f-iEElogtanhf^i^-^) (3.103) 

h=l h=l h=l k=l 



+ F3^t 

J — OO i 1 



l sinh(y[ 1) - y) 

being the integral evaluated on a zero we do not need to prescribe the principal value, according 
to (|3.97|) . A similar manipulation in the other strips leads to 



h=l 



^-^(V^)) (3-104) 



\ E »f * " i E E A i-' E logtanhf?!^- - 1 J 

h=l h=l r=l k=l 



/oo j L—1 m 9 ^ 

-oo 1 l_i smh in; 



, i sinh(y[ 9) - y) 



q = 2,...,L-2 
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We add all these expressions for different strips q = 1, . . . , L — 2 

mi L-2 m q . L-2 rn q m r -(g) „(r) 

-Y,^ = ^ E E - i E 4w E E lo S tanh(^^ - i£) (3.105) 

h=l 9=1 h=l 9,r=l fc=l fe=l 



We can use the following identity to evaluate the second term on the right hand side 

7T 7T 

logtanh(a; — i— ) +logtanh(— x — i— ) = —in (3.106) 
This identity is true in the fundamental determination of the logarithm and yields 

mi 1 L-2 m 9 £_3 



^E^ = ^EE^E-w ( 3 - 107 ) 

fe=l g=l h=l 9=1 

[°° dy I S , r , , S 



, 9,r=l 



h=i sinh(y^ - y) 



In the expression for the energy (|3.82j) we recognise the sum of exponentials computed in (|3. 107|) . 
We substitute it and we obtain an algebraic term and an integration term: 

-. L-2 ra q L-3 g 

E = ( 3 - 108 ) 

9=1 h=l 9=1 

The integral term is defined as 

5 = 2/ dy I Ae-yL^y) - £ A q>r L r (y) £ — } (3.109) 

J-qo I \ l x sinh(y^ — y) / I 



9,r=l 



and it will be evaluated later with Rogers dilogarithms. The first and second of the algebraic terms 
can be easily worked out using the relation (|3.74p between the two families of quantum numbers: 



L-2 m q L-2 m«j _ L-2 



^EE4 9) = EE^4E^- ( 3 - 110 ) 



2 / j / j ft / j / j ft 2 

9=1 h=l 9=1 h=l 9=1 



Using the Cartan matrix (|2.ip the following identity is easily proved: 

L-2 , L-3 



o E m l ~ o E m 9 m 9+l = T mT ^ m (3.111) 



2 ^ <? 2 ^ H ^ 4 

9=1 g=l 

where we use the same notation as in section 12.31 Using (|3.110|) and (|3.11ip in (|3.1U8p we obtain 
the following expression where the quantum numbers and the content of zeros enter explicitly 

, L-2 m„ 

i^m^m + EE^-^- (3-112) 

9=1 h=\ 
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To evalutate the integral contribution S in (13.1090 . it is convenient to compute the following ex- 
pression in two alternative ways 



Sn 



dy 



■ log(l + e 9 (y)) - log — 



(3.113) 



dy <- - ■> dy 

The absolute value here is convenient because it avoids us keeping track of the imaginary part that 
could appear when e q (y) is negative; this imaginary part does not contribute to the final value so 
we safely remove it from now on. Moreover, the full integrand is always real because of the reality 
condition (|3.96|) . The quantity here defined will be evaluated in one way as sums and differencies 
of dilogarithms and in a second way by using the TBA equations. 



3.3.6 Evaluation with TBA equations 

We evaluate the whole sum 

L-2 

E 5 « ( 3 - 114 ) 

q=l 

and we will see that many terms will cancel. We need to insert in (I3.113P the function log |e 5 (y)| 
and its derivative as obtained from the TBA equations: 

dktWI - + E^(E —±-<z- + K.K<y)). (3.H5) 

=1 smh(y - yJr) / 



dy 



r=l \k= 

One can easily realize that, within the whole sum (|3.114|> . all the terms that contain a convolution 
mutually cancels. Indeed, each term of type 

d log\e r (y)\ 



dy 



L r (y) 



(3.116) 



contains a product K * (L' r _ 1 + L' r+1 )L r and each term of type 

-log\e r+1 (y)\ L' r+1 {y) (3.117) 

contains —K * {L r + L r +2)L' r , ^. In these triple products the convolution and the ordinary product 
can be exchanged and this produces the cancellation of all the convolution terms. 
From the terms like f|3. 117|> we also have contributions of the form 



tanh 



*(r±l) 

y-y k 



dy (-l)log 

i 

where the absolute value enters because of f|3. 113[) . Integrating by parts we obtain 



(3.118) 



(-l)log 



tanh 



y-y k 



dy 



dy (-1) 



sinh(y - y| r±1) ) 



L r {.y) 



(3.119) 



sinh(y^ r±1) - y) 



where the first contribution vanishes. The remaining term in (13.119j) sums with a analogous con- 
tribution from (I3.116P and we are left with the following expression 



L-2 
q=l 



dy 



L-2 



h=l smh(^ - y) 



,q,r=l 



s 



(3.120) 



that precisely matches the quantity (|3. 109|) that enters the energy. 



36 



3.3.7 Evaluation with dilogarithms 

We start the evaluation of the first line in (|3.f 13|) by performing a change of the integration variable 



t = e«(y)>-l, dt = dy[e«(y)}' 
so that S q is uniquely fixed by the asymptotic values of e q 



S„ 



£9(+00) dt( l ° g{1 + t) log|t| 



Ui(-oo) ' t 1+t 

The integration interval can be modified passing through the point zero 



(3.121) 



(3.122) 



S a 



£g(+oo) dt f log(l + t) log |*| 
V t 1 + t 



gg( -° o) / log(l+t) log |*1 
u V t 1+t 



(3.123) 



so we need to evaluate single constituent blocks given by the integration from zero to a certain 
asymptotic value t a > —1. We introduce the function L + to label the main block 

+ y ' 2 J V t 1 + t 



1 f 1+a ^ dy / log(l - y) + logy 



o 



y 



c 



if t a ^0, 



1 

2./,, 



~* a ^ / 1 °g( 1 ~ y) + lo sy 



i + V' 

C(-t a ), if ia^O 



(3.124) 



where we performed simple changes of the integration variable and used the dilogarithm function C 
defined in Appendix A. Notice that, so far, we only used the dilogarithm function in the fundamental 
interval [0,1]. 

The integral under investigation (13.1130 can finally be expressed as the difference of the con- 
stituent blocks in (|3.123p through the function L + 



S q = 2 



L + (6*(+oo))-L + (e*(-oo)) 



(3.125) 



3.3.8 Excitation energies 



Our purpose is to make explicit the expression (|3.125|) . Using the known identity ()A.5|) in the case 
t ^ we obtain 

t 



L+(t) 



1 + t 



6 \l + t 



Similarly, using the continuation of the function (|A,2j) we have 



L+(t) = -C(-t) 



7T 



1 



6 '~U + t/' 

so we can summarize these two cases into a single expression 



t > 0. 



t ^ 



1 

2./,, 



t 



i+t 



6 \l + t 



(3.126) 



(3.127) 



(3.128) 
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Using now the asymptotic values (|3.20p . (13.930 . the argument of the dilogarithm in (13. 128j) takes 
the quadratic form 

1 sin 2 r sin 2 9 

= — o — or — — 3.129 

l + t simV sin 2 (g + l)0 

respectively for the e q (— oo) or e 9 (+oo) case. Finally, using (|3.128|) . the sum of expressions (|3.125|) 



can be rewritten as 

V-4e(m*(+~))-m-(-~») -4e[<^) -4 ■ f°,J ] < 3 - i3 °) 

7r z ir z ' V / 7r z ' L Vsnrar/ Vsim(o + l)#/J 



where we introduced the factor 3/7r 2 to make a comparison with the dilogarithm identity stated in 
(JXT]). Indeed, the tp in the first term specializes to ip = r = J so that 



r , -2 . r L-2 . 2 

o ^ — v _ / sin r \ o v— v _ / sin r 



U V — V _ / Sill 7 \ U V — V _ / Sill 7 \ , . 

8{p,2,L-2) = -=yc(- r - T , 7— )=^Y£[— -2—) 3.131 

v ' 7r 2 ^ Vsin 2 (g + l)r/ vr 2 ^ Vsin 2 gr/ 



g =i ^ ' 9 =i 

and the last equality is justified because 

,2^ \ ^2 



„ / s „ / sin r \ -/r . 
£ 1 = £ , = -. 3.132 

Vsm (L - 1)tJ o 



In the third term of (|A.7jl we have <p = 9 = so that 



6 4-^ „ / sin 2 9 \ 6 v— ? „ / sin 2 \ , 

s(0, 2, L - 1 = V C( — — = — V £ — =- — + 1 3.133 

V ; 7T 2 ^ V s in 2 ( 9 + 1)0/ vr 2 ^ V s in 2 (g + 1)0/ V ' 



where we used 



„ . , / sin 2 9 \ ir 2 , 
£1 =£(_-) = —. 3.134 
Vsm Lt// o 



With the results (|3.13ip . (j3.133H . the expression (j3.130n becomes 

3 6 

-^S = s(0,2,L-2) + l-s(0,2,L-l) = c L = 1 - — — — (3.135) 

7T Z + lj 

where the identities (|A.7|) . (|A.8|) were used and the vacuum boundary conditions were specified 
r = s = 1. We can insert this result in the expression for the energy (|3.112|) to obtain 

L-2 m q 

£ = -g + ^ r Cm + ££4 9) (3.136) 

q=l k=l 

and we see the appearance of the central charge of unitary minimal models. We have obtained a 
completely explicit expression for the energy of the states of the theory. 

The identities (|A.7|) , (|A.8|) remind us the analogous expressions for the conformal dimensions of 
fields in unitary minimal models. Indeed, the full TBA calculation of this sector can be immediately 
extended to all cases. More difficult is the task of deriving the (m, ra)-system that applies to an 
arbitrary sector (r, s). We do not discuss this issue in the present paper. 
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4 s£(2) 1 WZW Lattice Models and XXX Hamiltonian 
4.1 Definitions 

The bulk face and boundary triangle Boltzmann weights of the critical level-1 WZW vertex models 
in the vacuum sector are 



l-z. 



<-K- 



1 



(4.1) 



(4.2) 



where z = e lu is now an multiplicative spectral parameter. Commuting double row transfer matrices 
can now be defined following [35]. The Hamiltonian limit of this model gives the XXX quantum 
spin chain. The quasiparticle description of this model, including the (m, n) system, tower particles, 
string patterns and Cartan matrices, is given in Section 2. In the rest of this section we consider 
the functional equations and solution of the resulting TBA equations. 

4.2 Bethe ansatz 

The Bethe ansatz equations satisfied by the double row transfer matrices of the XXX model are 
given by 

zD(z)Q{z) = (z + ±) 2N+1 Q(z-l) + (z-±) 2N+1 Q(z + l) (4.3) 

where 



Q(z) = H(z- Zj ) 

3=1 

and the physical analyticity strip of T(z) is —1/2 < Hez < 3/2. 



(4.4) 



4.3 Functional equations 

The functional equations given in section [3~2l hold also in the present case, with the difference that 
here there is no truncation so the fusion level q takes all integer values up to oo. Moreover, the 
parameter A is now an arbitrary real number. We are interested in the limit A — > so it is useful 
to rescale the spectral parameter and take the limit 



u 



A 



uX 

IT 





(4.5) 
(4.6) 



This rescaling and limit transform all trigonometric functions into polynomial functions of the 
rescaled variable u. We introduce D (u) = Dq and 



D\ = D g (u + kir), Q k = Q(u + kir), s k {u) = 2- + k, f k = (-l) N (— + k) 



2N 



(4.7) 
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The definition (13.101) becomes 

dl= q , f ? . K?. (4-8) 

S-2S2qJ-lfq 

and the Y-system has a lower closure only 

dgdf = (l + d^Hl + d?- 1 ) (4.9) 

= (4.10) 

We also rewrite (I3.9P for g = 1 as 

S - ,Sl -DlD\ = I + dl (4.11) 



S-2S2qf-lfq 

We remark that the periodicity of it that was a feature of all the transfer matrices (|3.13p is now 
lost after the rescaling of the spectral parameter. The new analyticity strips for D q (u) and d q (u) 
at each fusion level in their respective analyticity strips are 

- 1 7r < Ren < — 7r, - g ^ 1 7r<Ren< 3 2 ^ ^ ( 4 - 12 ) 
Defining the shifted transfer matrices 

D q (u) =£> 9 («+^7r), d 9 (u) = d"(n + ^vr) (4.13) 
it follows that these transfer matrices have the common analyticity strip 

- 7r < Re(u) < 7r (4.14) 
and satisfy the same crossing symmetries 

£>«(«) = £> 9 (-«), = d 9 (-u) (4.15) 
In terms of shifted transfer matrices, the TBA functional equations take the form 

d\u - \) d\u + I) = (1 + d q -\u)) (1 + d q+ \u)) (4.16) 
We evaluate the asymptotics taking the limit A — > of (|3.2Up and we get 

d q (+ioo)=q(q + 2) (4.17) 
Observe that now the asymptotics are the same for all sectors. 

4.4 Analyticity and solution of TBA 

Most parts of the derivation of the TBA equations for the truncated case carry on here, without 
modifications, so we limit ourselves to the results of the calculation. 

Starting with the energy calculations, the bulk and boundary factors survive in the limit (|4,5p 
so we can remove order N and order 1 zeros from D 1 as in (|3.24[) 

dLM = t — r^m — = ^Lte(«-£). «e(~,§ir). (4.is) 
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We easily get 

£finite(«-f) Dl nite (u+^) = l + d l (u), U€ (_I,|). (4.19) 

We now remove the zeros of the level=l transfer matrix, that we denote as 

DU e (±w£ ] ) = D l {MyV) = , y« > . (4.20) 
To that purpose, we construct the following function 



/c=i 



(i) , ■ (i) • ■ 

tanh^ tanh^ 



uG(-7r,7r) (4.21) 



that guarantees the complete removal of all zeros from the analyticity strip. Again, we rotate the 
variable 

e q (x) = d q (ix) (4.22) 



and use the analytic and non-zero function 

Amite («) 



A{x) 



Z(u) 

to write the functional equation 



(4.23) 



A(x + i'-) A{x-i'-) = l + e\x), UG (-|,|). (4.24) 

We solve it by Fourier transform, as in (|3.36p . and, at the isotropic point, the same equation 
obtained in (|3.39[) . 



2t , , , . 1 

— E + higher order corrections m — 



It is now clear that we can directly rescale and take the limit (|4.5p in the equations for the truncated 
case and we will get the equations for the untruncated one. From (|3.52|) we have 

loge q (x) = 8 x , q N log tanh 2 | + log tanh 2 | + C {q) (4.26) 

oo m j (j) (J) oo 

+ log n n [ tanh x -^- tanh x -^-\ Aq,j + e a. j ( k * mi + ^)) w 

i=l fc=l 3=1 

with the same kernel as in (I3.37p . We evaluate the +oo limit of the equation to fix the integration 
constant. Using (13.530 we get 

1 oo 

loge^(+oo) = - J2 A qJ log(l + e^(+oo)) + C® (4.27) 
i=i 

and the asymptotic values (|4.17p force the vanishing of all integration constants 

C^=0. (4.28) 
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The zeros are fixed by the untruncated family of $ functions 



7T X IT X IT 

■$i(x) = i\ogs q {x - i— ) = iSi t q A r logtanh 2 (- - i— ) + ilogtanh 2 (- - i— ) + 



+ i ^ A Q,j [ lo § tanh ( 



3=1 k=l 

OO 



(i) (i) 
^-5- - t-) + log tanh( ^-S- - »- 



00 Iog(l + 

-oo 



27T sinh(x — y) 

with the quantum numbers 

*9( y («)) = 7T n < 9 > 4 g) = 1 + 2(4 5) +m q - k) . 
The scaling limit is computed precisely as in section 13.3.41 The result for the energy is 

E = -—\ ^ e -^log(l + e 1 (y)) + 2V-e- 
with the following set of TBA equations 



k=l 



1 



log 



-4 <J 1)(J e x + lo; 



oo r»j 

:nn[ 

3=1 k=l 



tanh 



+ 



+ ^4 Z)i (#*log(l + ^'))( 



9 = 1,... 



and the following quantisation conditions and auxiliary equations to fix the zeros 

oo m r ~(r) _ 

&(x) = 45 llff e-* + tX^ ir X;iogtanh(^=^--4) 



r=l fc=l 



9i r 



(?) 

7T 7l£" = 7T 



dy 



\og(l+e r {y)) 
2ir sinh(x — y) 

1 + 2(1$ + m q — k) 



along with the constraint on particle composition given by 

mi m2 ^ • • • ^ rrtfc ^ • • • ^ 



(4.29) 



(4.30) 



(4.31) 



(4.32) 



(4.33) 



(4.34) 



After, the — oo limit can be evaluated from the TBA equations or simply by taking the limit of 
(PEMD for A -> 



-oo) 



1 



(4.35) 



as in (I4.17P we see that the asymptotics do not distinguish the sector. 
4.5 Exact energy from excited TBA 

Our goal is to get an explicit expression for energy eigenvalues. The method is very similar to the 
one introduced in section [3.3.51 but we need to introduce slight modifications to avoid intermediate 
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divergent sums as those appearing in (j3.107H . Moreover, the non-cancellation of a number of terms 
compared to the truncated case requires to introduce a very compact notation to keep expressions 
within a reasonable length. We define 

-(<?) _ -W 

= ~^logtanh( ^ 2 Vk (4.36) 
and observe that the identity (|3.1U6|) takes now the form 

Bqh,rk + Brk,qh = ~ „ (4.37) 

Also, we define 

L r (y) 

47r 2 sinh(yj^ — y) 



G r , qh = / dy — '^f- (4.38) 

Th 

here the integrand is always finite because the numerator vanishes on the adjacent strip zeros as in 
(I3.97P therefore the principal value is not necessary. We rewrite the quantization condition (14.330 
for a zero in a generic strip y^' as 

h^e-*™ = + Kr^lKrh + E A i,rG r , qh (4.39) 

r,k r 

For both B and G notations, the comma separates variables that refere to different objects, for 
example in G r ^ q h the label r referes to the numerator while qh to the denominator. In the sums, 
an index as in ^ r means to sum on all possible values of r; which are its possible values it is clear 
from the context, for example the following sum 

oo m r 

E B Qh,rk = E E B 1 h ' rk ( 440 ) 
r,k r=l k=l 

has no other possible interpretations because r ranges on positive integers and k, according to 
(|4.36p . labels a zero so it can take values in the corresponding strip r only. Also, no sum on 
repeated indices will be used. 

In (|4.39|) we sum on all zeros of all strips up to a given Q > 1, namely on all q = 1, . . . , Q and 
on all h = 1, ... , m q . Actually observe that the contribution of a strip q > 1 to the right hand side 
is zero so that sum is independent of Q. The sum is reorganised 

n n n 

2 ^ 



= \i2 n h ] + E A c,rB qh ,rk+ E A q>r G r , qh = (4.41) 

q=X,h q=l,h,r,k q=l,h,r 

- Q Q Q-l 

o E n h + E B <lh,q-lk + E B qh,q+ik + E B Q h ,Q+ lk + 



q=l,h q=2,h,k q=l,h,k h,k 

Q-l Q 

+ E Gq+l,qh + E Gq-l,qh + E Gq+1,QH 
q=l,h q=2,h h 

In the second line, the sum on the quantum numbers is evaluated with (|4.33|) and gives 

Q Q 1 Q 



q=l,h q=l,h q=l 
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the second and third terms are evaluated with (14.370 and yield 

2 



1 Q-l 

- m q m q+1 (4.43) 



=1 



The sum evaluated in (|4.4ip can now be inserted in the energy expression (|4.3ip so we write 
Q 2 i Q- 1 

E = ^2l { h q) + ^ + -Ym q+1 (m q+1 -m q )+ (4.44) 

q=l,h q=l 

Q O-i ,oo dy 

+ E &Qh,Q+lk + E G q-l,gh+ E Gq+lah + Y^GQ+l&h- \ V L l(y) = 

U 1. _ r, i _ -I . r_ J — OO 7I_ 



q=2,h 9=1, ^ 

r(9) 



5=1, h g=l 



and the newly introduced quantities are defined as 

S Q = 8 dye-y^-S^Y.G.-i^-S^Y.G.+^h (4.45) 

g=2,/i g=l,/i 

T Q = E^ ( 4 - 46 ) 

We need to sum up to Q the values 

/OO 
-Re(e q (y)yL q (y) + Re(e q (y))L' q (y) (4.47) 
-oo L - 1 

where the real part is taken to remove all imaginary contributions coming from the zero terms in 
(P~32|) . namely 

_ -M _ 
Relogtanh ^ ^ k ) = log | tanh( ^ ^ k )| (4.48) 

With some integrations by parts in (|4.47j) we get the following expression 

/OO 
dy e- y Lx (y) - 8vr 2 A 5 , r G q , rk + E 4,,r {L q ■ K * L' r - L' q ■ K * L r ) (4.49) 

where the dot product • indicates the bilinear form 

/OO 
dyf(y)g{y) (4.50) 
-OO 

and the convolution "*" has higher priority than the bilinear form Actually, the two types of 
product can be exchanged because the kernel is even, for example 

L'-K*L r = L r -K*L'=L'*K- L r (4.51) 
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This property is extremely useful in the following sum because all bilinear terms cancel except 
those with largest indexes 



3=1 



/oo ^ 
dy e- y L x (y) - Sir 2 ^ G q ^. qk - 8vr 2 ^ G g+1 , qk + 

"°° q=2,k q=l,k 

- 8tt 2 ^ G Q,Q+i k+K * L' Q+1 ■ Lq — K * L Q+1 ■ L' Q = 
k 

= S Q - 8tt 2 G Q:Q+1 k + K* L' Q+l ■ Lq- K * Lq+i • L' Q 



(4.52) 



and we see that the quantity appears in this sum; in this way, the sum of (I4.47P enters into the 
energy expression. The evaluation of (I4.47|) with dilogarithms has been done in section 13.3.71 and 
holds in the present case; in particular we will use (I3.125p . From the asymptotic values (I4.35|) and 
(|4.17p . using also the properties given in Appendix A, we get 



L + (e q (+oo)) = C 



£ 9 (+Oo) 

1 + ei{+oo) 



C 



g(q + 2) 
(q + l) 2 



L + (e q (-oo)) =£ 



and finally 



S q = 2 



g"(-oo) 
1 + ei(-oo) 

1 



C 



q 2 - 1 



6 V(g + 1) 2 

2 



c 



c 



(9+1) 



In the sum on q all terms disappear but two 
" " 1 



q=l q=l 

and (I4.52|) reads 



Eft? 



(9+1) 



2£(1) - 2C 



(Q + i) 



r( 1 

3 " 2£ l(Q + i)2 



(4.53) 
(4.54) 

(4.55) 
(4.56) 



<?Q = 8vr 2 ^ - K * L' Q+1 • L Q + K * • L> Q + ^ - 2£(^ T 



(4.57) 



Now we are ready to express the energy. From (I4.44p . using also H4.52H and ()4.56p . we have 
E = -^ d +E^ + i- + lj:rn q+1 (rn q+1 -rn q ) + UQ 



q=l,h 



q=l 



1 



C 



1 



8tt 2 24 4tt 2 V(Q + 1)2 



(4.58) 



+ ^2 ( K * L Q+l - L Q~ K * L Q+l ■ L 'q) + BQh,Q+l 



h,k 



By construction, the energy expression H4.58H is actually independent of Q because in (I4.39P the 
contributions with q > 1 actually vanish. We will show that in the limit Q — > oo the energy (I4.58|) 
will take a completely explicit form; to argue this, we need to know the behaviour of the transfer 
matrix eigenvalues for very large strip index q. 
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The number of 2-strings n\ grows thermodynamically as N while all the other one- and 2- 
strings numbers remain finite therefore the first line a = 1 of (|2.10p becomes the trivial identity 
00 = 00. This mechanism applies to all states; suppose a state is described by a lattice path 
of lenght No and maximal height Lq = [Nq/2\+1. As we increase N, the same path fits in the 
extended lattice with just the addition of particles of type n\ (Dirac sea) on the right-most part of 
the lattice. All the other particle numbers are unchanged. Explicitly, the following particle contents 
do not vary with N 

m a for a = 1, 2, . . . , Lq — 1 

n a for a = 2, 3, . . . , L - 1 (4.60) 
n 

while the newly- arrived particles m a , n a , a = Lq, Lq + 1, . . . are fixed by consistency as 

m a = no for a = Lo,... (4.61) 
n a = for a = Lq, . . . 

where we do not specify an upper limit because these equations hold true for all the following strips. 
It is also useful to define 

rrioo = lim uil = no (4-62) 

L^oo 

If we suppose that (|2.10p is verified for a size No = 2(Lo — 1), the highest tadpole equation a = Lq — 1 
becomes (remember that mi r i = uq = rriL ) 

m Lo - 2 + mLp-i m Lo _ 2 + m Lo 
m Lo _i + n Lo _i = = (4.63) 

therefore the last member takes the form of an A-like adjacency rule. With the given definitions 
(|4.6ip . for all L ^ Lq the following A-like equality holds 

m L -i+m L+1 

m L + n L = = n (4.64) 

We unify all this in the (m, n)-system notation 

00 

m a + n a = \ ^ A a>b m b a = 2,3,... (4.65) 

6=1 

where the adjacency matrix is now an infinite ^-like matrix. In summary, the identity rriL —i = 
no = mi for all L ^ Lq embeds the tadpole diagram in an infinite diagram. The property 
(|4.6ip is extremely important: it states that high-numbered strips have the same particle content, 
from a given Lq on. The rectangle Ao = 2(Lo — 1) times Lq represent the smallest lattice at which 
a generic state (|4.60p exists. 

With (|4.6ip we can safely take the limit Q — > 00 in the energy (|4.58|> . In that limit, the sum 
on quantum numbers is finite because strips (labelled by a) do not contain 2-strings for large a; 
the sum on the m-family particles is also finite because m q+ i — m q = for all q ^ Lq — 1. The 
dilogarithm vanishes because its argument vanishes. 

We are left with the following sum of terms to be evaluated in the limit Q — > 00 

( " Yl G Q,Q+i k + G Q+i,Qk) L' q+1 -L q -K* L q+1 •L' q )+^2 B Qh>Q+1 k (4.66) 



h.k 
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We assume that the full information content of the high strips is the same namely that the following 
limits exist at large q 



V (<?) (°°) 

lim e q (x) = e°°(x) 

9— ►oo 



(4.67) 



The two conditions are not fully independent because of the TBA and auxiliary equations. This 
limit says that the high strips coincide, they carry the same information (see the section on numer- 
ical evaluations for a justification). With this assumption, we can evaluate (|4.66p by substituting 
the full expression f|4.38|> 



lim 

Q^oo 



>+l,Qk 



Q— >oo 



-L Q (y) 



+ 



U 



i(y) 



sinh(y[, Q+1) - y) sinh(y^ - y) 



The convolution term is evaluated by remembering the commutativity (14.511) 



(4.68) 



(fee 8^2 { K * L ' Q+1 ' L Q~ K * L Q+i " L 'o) = ( K * L 'oo ■ L °c - K * L oo ■ L' c 



(4.69) 



The last term must be written explicitly and, according to (|4.62p and to (I4.6ip . we will use iuq 
nice because Q is sufficiently large 



• moo ~(Q) _ -(Q+i) 

lim Y, B Qh,Q+±k = lim ~ £ logtanh(- % ' 

Q^oo * — ' Q— >oo Z7T ' — ' 



h,k 



h,k=l 



nice ^(oo) _ ~(oo) 

' ^logtanh( ^ y x -4; 



2vr 



h,k=l 
rrioo 



• 7^ ^ 
*4 } 



2tt 



log tanh(— i— ) + 



/i=l 



- y 

2vr ^ 



h,k=l 
h<k 



~(oo) ~(oo) _ ~(oo) ~(oo) 

logtanh( ^ - il) + logtanh( ^ - i£) 



moo 1 n x 1 2 

-— - ^ooK - 1) = -^ m oo 



(4.70) 



where we have used the identity (|3.106p . Notice that we always use the logarithms in their funda- 
mental determination. In (I4.58P all terms have been evaluated in the limit Q — > oo so we write the 
final result 



oo rn q 



9=1 h=l 

oo rn q 



E = + E E 'f + y + ^ E "it 1 '"!*! -'».)- j m 



9=1 h=l 



9=1 



2 

4 --oo 



(4.71) 



where C is now the infinite A^-like Cartan matrix, in agreement with the (m, n) system (|4.65p . 
We have to compare the energy expression (|4.71|) with the finitized character expression (I2.3U|) . As 
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noticed in (|4.60p . every state is characterized by Lq namely the minimal lattice that contains it. 
Considering a finite lattice and matrices and vectors of length Lq, it is easy to prove that 

^m T C t&d m = ^rn T Crn-^ml (4.72) 

where C ta( j is the tadpole T 1 ^ -y Cartan matrix and C is the Al _\ Cartan matrix. Observe that, 
a priori, on the left hand side the limit Lq — > oo is not defined because we don't know how to 
treat the entry corresponding to the loop in the diagram, namely the entry (Lq,Lq) of the Cartan 
matrix, while in the right hand side the limit is perfectly defined thanks to the properties (|4.6ip . 
More explicitly, in the right hand side the general term takes the form 

m q (m q -m g _i) ^ ^ 

and it is trivially zero from Lq on because of (|4.61[) . 
4.6 Tadpole truncation for numerical TBA 

In (I4.67P we have formulated an assumption on the behaviour of the transfer matrix eigenvalues 
in strip q, with q growing to infinity. Here we try to motivate this assumption on the basis of 
numerical solutions of the infinite TBA system of (|4.32p , (j4.33j) . 

It is important to remember that a strip is directly related to the neighboring strips only, with 
an adjacency rule, therefore we face the problem of choosing a suitable truncation. Inspired by 
the argument given near (|4,65p and also by the assumption in (|4.67p . we use a tadpole truncation, 
in which the amount of information that should come in the last strip Q from the strip Q + 1 is 
replaced by the strip content of Q itself (this Q is not the same of Section |4.5|) . In other words, 
(|4.32|) . (|4.33|) are unchanged for all strips q < Q but for the last strip we replace them with 

Q m j _ r (i) Q 

\oge Q {x) = log Yl H tanh X ^ + ^ (K *\og(l + e j ))(x) (4.74) 
j=Q-lk=l ~ j=Q-i 

r=Q-lk=l r=Q-l J 00 K y> 

Quantization conditions and (|4.31|) are unchanged. This truncation can be equivalently formulated 
by saying that the adjacency matrix of in (j4.32H . ()4.33p is replaced by that of the tadpole Tq 
as in Section [2] and expecially in (I2,10p , Notice that the counting in (I3.10ip makes inconsistent a 
naive truncation in which the last strip Q contains the previous one Q — 1 but not the Q itself: 
that would make no space to allocate the required number of 1-strings. 

We have implemented numerical TBA equations for the truncations T' 2 , Tg, and the actual 
numerical results fully confirm our assumption. For the second state in Figure [10} characterized by 
two zeros in each strip, with quantum numbers 

lW=/W=0 V q, (4.76) 

we obtain the zeros positions in Table HJ The distance between corresponding zeros in different 
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strips \yf' — yi | can be evaluated from Table [Hand suggests a vanishing behaviour 





i = 2 


3 = 


= 3 


J = 4 


\y? 




0.6038 
0.3949 


0.3S 
0.2* 


361 

m 


0.2890 
0.2323 



consistent with the existence of the limit (|4.67p . The actual positions grow with j but we expect 
this growth will slow down later on. From Table [T] we also observe that the energy value decreases 
toward 1, as expected for this state, see Figure [TU1 

The result in (|4.77|) is precisely as expected: different strips tend to resemble each other, with 
the same zeros and the same eigenvalues (plots of the transfer matrix eigenvalues are not given 
here). We obtained similar results for various other states, for example for the first of Figure [TU 

5 Discussion 

In this paper we have introduced a combinatorial formalism, based on paths and quasiparticles 
(particles and dual-particles), to classify the eigenvalues and eigenstates of the transfer matrices of 
critical two-dimensional Yang-Baxter integrable lattice models. This is achieved by implementing 
a conjectured energy-preserving bijection between patterns of zeros (in the plane of the complex 
spectral parameter u) of the transfer matrix eigenvalues and lattice paths coinciding with one- 
dimensional configurational sums. The TBA equations are solved for the conformal spectra of 
energies using the analyticity information encoded in the patterns of zeros, that is, the relative 
locations of the 1-strings (dual-particles) and 2-strings (particles). This program has been carried 
to completion for the ABF RSOS and XXX models in their vacuum sectors. However, we expect 
our methods to generalize to other sectors and to other models including the D and E RSOS models 
as well as parafermions. Finally, our methods should extend to boundary conditions other than 
fixed boundary conditions such as periodic and toroidal boundaries. 
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Table 1: Zeros and energies from the numerical solution of the TBA system at three different levels 
of the tadpole truncation are given here. 





E 


y? 




(2) 

y\ 
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(4) 

y\ 
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y 2 


rp/ 


1.00714 


-0.2783 


1.6430 


-0.6732 


2.2468 










rp, 


1.00397 


-0.2874 


1.7675 


-0.8257 


2.6191 


-1.1155 


3.0052 
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1.00253 


-0.2911 


1.8293 


-0.8989 


2.8152 


-1.3431 


3.4124 


-1.5754 


3.7014 
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A. Dilogarithm Identities 

The Rogers dilogarithms are defined as 



^)=4r^f^ l +^ 1 o<x<i. (a.i) 



2 7 v y 1 - y 

Following |43[ [33] , we extend the definition to the whole real axis by 



7T 2 

C(x) = —-Cix^ 1 ), if x>l, (A.2) 

C{x) = £ (tz^) ~Y> if x < °> (A ' 3) 

£(0) = 0, £(1) = ^, £(+oc) = ^, £(-00) = -^. (A.4) 

DO 6 

We will use the following identity that remains true for all x 

C{x) + C{l-x) = ^r (A.5) 
6 

In [33], an important identity is stated and proven. We write it with special reference to our case. 
Following eq. (2.1) of Kirillov, we define 

s(j,2,fc) = — > £( . 2 j, y= . (A.6) 

Corollary 2.5 of Kirillov states that 

s(r - 1, 2, L - 2) + 1 - s{s - 1, 2, L - 1) = c L - 24/i r . s + 6(r - s)(r - s + 1) (A.7) 

where 

6 , \(L + l)r - Ls} 2 - 1 , K , 

and L, r, s are certain positive integers. 
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